Metal-Insulator Transition in Two- and Three-Dimensional Logarithmic Plasmas 
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We consider the scaling of the mean square dipole moment in a plasma with logarithmic in- 
teractions in a two- and three-dimensional system. In both cases, we establish the existence of 
a low-temperature regime where the mean square dipole moment does not scale with system size 
and a high-temperature regime where it does scale with system size. Thus, there is a nonanalytic 
change in the polarizability of the system as a function of temperature, and hence a metal-insulator 
transition in both cases. The relevance of this transition in three dimensions to quantum phase 
transitions in 2 + 1-dimensional systems is briefly discussed. 

PACS numbers: 64.60.Cn,71.30.+h,52.65.Pp 
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Gauge theories where matter fields are coupled to com- 
pact gauge-fields in 2 + 1 dimensions have recently been 
much studied as effective theories for Mott insulators 
with competing orders HSU. Compact 17(1) gauge the- 
ories in 2+1 dimensions are capable of sustaining topolog- 
ical defects in the gauge fields in the form of space-time 
instantons Q. Such theories generically offer the possi- 
bility of featuring confinement-deconfinement transitions 
associated with a proliferation of topological defects in 
the gauge sector. It is hoped that phase transitions in 
matter coupled such theories 0, E| may be connected to 
difficult problems in condensed matter physics, such as 
the breakdown of Landau Fermi liquid theory and pos- 
sibly also spin-charge separation in strongly correlated 
systems at zero temperature in two spatial dimensions 
|2j. One such model considered recently is the compact 
abelian Higgs model where the gauge field is coupled to 
matter fields with the fundamental charge |8(- By itself, 
the compact U(l) gauge sector may be mapped onto the 
2 + 1-dimensional Coulomb gas with 1 /R- interactions be- 
tween "point charges" i.e., the instantonic topological 
defects of the theory As shown by Polyakov, the lat- 
ter system is always in a metallic phase (when the instan- 
tons are regarded as electric charges) More recently, 
it was shown that when matter carrying the fundamen- 
tal charge was coupled to the compact gauge sector, and 
critical matter field fluctuations were integrated out, the 
system could be mapped onto a gas of point charges in- 
teracting with a potential — ln(i?) in 2 + 1 space-time di- 
mensions, instead of 1/R Q- Using RG arguments, it was 
demonstrated that the system with the 2 + 1-dimensional 
logarithmic interaction may undergo a finite-temperature 
phase transition driven by the unbinding of dipole con- 
figurations 8], from a low-temperature dielectric regime 
to a high-temperature metallic regime(again when the 
instantons are regarded as electric charges). Here, we 
demonstrate this using Monte Carlo (MC) simulations. 

Whether or not there exists a low-temperature dielec- 
tric regime separated from a high-temperature metallic 
regime in a three-dimensional system of point charges 



interacting with logarithmic interactions and overall 
charge-neutrality, is presently a matter of debate in both 
the condensed matter and lattice gauge theory literature 
H 0- E3 • The theorem that dipoles cannot screen the 
Coulomb potential (in any dimension) 11] does not ap- 
ply to the case of a logarithmic pair potential between 
point charges in three dimensions. It is due to this in- 
ability of the dipoles to screen a Coulomb potential that 
the Kosterlitz-Thouless (KT) transition is possible: the 
low temperature dielectric phase is always critical below 
the KT transition temperature. In three dimensions, on 
the other hand, Debye-Hiickel theory is essentially ex- 
act. Thus, starting from a system of charges, screening 
is such that the screening length cannot become infinite, 
and the system is always in the metallic phase. On the 
other hand, if one starts from a three-dimensional dipole 
system, the system always stays in the dielectric phase, 
since dipoles cannot screen. Once more, there is no phase 
transition. If some departure from the Coulomb poten- 
tial occurs in three dimensions, dipoles could conceivably 
be able to screen. A 3D logarithmic interaction could in 
principle be screened to a potential that decays with sep- 
aration between the charges |9(. If so, no metal-insulator 
transition occurs at any finite temperature, analogous to 
the situation in a three-dimensional Coulomb plasma 0] . 
If one thinks of the gas of point charges as a gas of in- 
stantonic defects in a compact gauge field, the lack of a 
metal-insulator transition due to non-standard screening 
by dipoles would correspond to permanent confinement 
in the compact gauge theory. 

We therefore consider a logarithmic plasma in two and 
three dimensions on the lattice using MC simulations, 
and study the polarizability of the system as a function 
of temperature. The screening properties of these sys- 
tems determine whether they are insulators or metals 
and are governed by the dielectric constant e, which in a 
low-density approximation is given by the polarizability 
p of the system, e = l+n^dP- Here is the dipole den- 
sity, and Qd is the solid angle in a d-dimcnsional system. 
Since the polarizability of the system is proportional to 
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the mean square separation (s 2 



I 2 ) between 



the charges constituting the dipoles, it is natural to fo- 
cus on (s 2 ) in order to investigate whether the system 
is a dielectric or a metal. We find a low-temperature 
regime where dipoles are tightly bound, separated from a 
high-temperature regime where they are unbound. This 
implies the existence of a low-temperature insulating 
regime and a high-temperature metallic regime, sepa- 
rated by a genuine phase transition. For comparison and 
benchmark purposes, we compute the same quantities for 
the two-dimensional logarithmic plasma, where a metal- 
insulator transition in the form of a KT phase transition 
is known to exist 0, 0] . 
We consider the Hamiltonian 
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where the sum is over all pairs of sites of a two- or three- 
dimensional periodic lattice and qi is the charge at site i, 
N = L d is the total number of sites, and a is the lattice 
spacing. We have subtracted V(f = 0) from the po- 
tential, since only neutral configurations will contribute 
to the partition function, and we work at zero chemi- 
cal potential for the instantons 0|. This has the usual 
Coulomb form for d = 2, but differs from it in the d = 3 
case: instead of the usual 1/R potential of the three- 
dimensional Coulomb gas (3DCG), the power 3/2 results 
in a logarithmic potential. Taking the continuum limit 
of Eq. J3J) for d = 3, we find 
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where Ko(x) is a Bessel function, A — 2n/a, and the 
Debye-Hiickel screening length m _1 was introduced as 
an infrared regulator. AsrrnO we have 



V(r) = —^[l- 7 -ln(Ar)} 



(4) 



where 7 is the Euler constant. Therefore, Eq. J3J) behaves 
exactly in the same way as the two-dimensional potential. 

Let us first turn to the 2DCG which we use as a bench- 
mark on our method of probing the metal-insulator tran- 
sition, before proceeding to the three-dimensional log- 
arithmic gas (3DLG). At low temperatures this model 
consists of tightly bound dipoles with the two opposite 
charges separated by a distance |rj — rj\ = s of the order 
the lattice constant. 



In a 2D system, it is well-known that dipoles begin 
to unbind at a critical temperature Tkt |13| and at 
high temperatures the 2DCG is a fully ionized metal- 
lic plasma, separated from the low-temperature dielec- 
tric insulating phase by a genuine KT phase-transition. 
Thus, we expect no finite-size scaling of (s 2 ) below the 
KT-transition, whereas we should expect (s 2 ) oc L a ' T ' 
with a(T) < 2 at higher temperatures. Using an intu- 
itive low density argument, neglecting screening effects 
[l3j |. we can calculate the behaviour of (s 2 ) to leading 
order in L, 
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L 2 



T <T KT 

T KT <T< 2T KT (5) 
2T KT < T. 



Hence, a(T) is zero for low temperatures and a mono- 
tonically increasing function of temperature just above 
Tkt- 

Including screening effects in 2D shows that this con- 
clusion still holds, however the temperature at which it 
occurs is determined by screening. This is ultimately re- 
lated to the fact that dipoles contribute only through 
a correction to the coefficient of the bare logarithmic 
pair-potential. Hence, in 2D the functional form of the 
renormalized potential is unaltered, only prefactors are 
changed. In three dimensions, it is far from obvious that 
dipoles do not have a much more disruptive effect on a 
bare logarithmic pair-potential. 

To adress the issue of scaling of (s 2 ), we will use MC 
simulations and finite-size scaling to study (s 2 ) both in 
the 2DCG as well as in the 3DLG. In our simulations, 
the particle number is not conserved. However, during 
the simulations the system is maintained electrically neu- 
tral. The MC moves involve creation and annihilation of 
charges, applying the Metropolis algorithm in this pro- 
cess. Starting at some randomly chosen lattice site, an 
attempt to insert a negative or positive charge at ran- 
dom at this site is made, with an opposite charge at a 
nearest-neighbour site. The move is accepted with prob- 
ability exp(-A£;/T) = exp[-(W new - H \ A )/T], and this 
is done in all d directions, before we move to the next 
site. It is clear that placing a charge on top of an op- 
posite one corresponds to the annihilation of the existing 
one. In order to measure s we have to keep track of which 
two charges belong to each other in a dipole, since there 
is no physical link defining the dipole. Two charges in- 
serted into the lattice at the same MC move is chosen as 
a dipole, and if one of these charges subsequently is an- 
nihilated, the effect of this annihilation is to diffuse one 
of the charges in the dipole. Hence, such a diffusion may 
increase or decrease s for this dipole |l5|. One sweep is 
defined as going through all the sites in the lattice once, 
and at every tenth sweep we sample s 2 averaged over the 
system. The sample is, however, rejected if there are no 
charges in the system at this MC time, since there is no 
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information on (s 2 ) in such a configuration. The ther- 
mal average obtained at the end of the simulation is thus 
taken only over non-empty configurations. 
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FIG. 1: Results from MC simulations of the 2DCG, where 
1.0- 10 5 sweeps were used at each temperature, a) (s 2 ) versus 
T for a selection of the simulated system sizes L =8, 12, 16, 
20, 24, 28, 32, 36, 40 and 48. Errorbaxs are smaller than the 
symbols used, b) a versus T found from fitting the data of 
(s 2 ) at different T to AL a where A is a constant. A selection 
of such fits are shown in c). We note that (s 2 ) is practically 
independent of L up to a certain T. 

The results from the simulations on the 2DCG are pre- 
sented in Fig. n A total of 10 5 MC sweeps at each tem- 
perature are used to produce the first plot. From the 



(s 2 )-data of system sizes L = 8, 12, 16, 20, 24, 28, 32, 36, 
40 and 48, we extract a for temperatures in the interval 
(1.0, 1.52), which is shown in the second plot. It is evi- 
dent that there are two distinct regimes of temperatures, 
one in which the charges of almost all dipoles are bound 
as tightly as possible, the separation of the charges cor- 
respond to the lattice constant. In the high-temperature 
regime the dipoles have started to separate, reflected by 
a scaling of (s 2 ) ~ L a ^ with the system size. The two 
regimes are necessarily separated by a phase transition, 
since in the low-temperature regime a(T) = while in 
the high-temperature regime a(T) ^ 0. This necessar- 
ily implies a non-analytic behavior of a(T). An attempt 
to determine the transition temperature from these plots 
yields approximately 1.32. This is slightly less than the 
early results of Saito and Miiller-Krumbhaar of 1.35 in 
our units [l6|. but is in excellent agreement with much 
more recent simulations by Olsson [ljj. It provides con- 
fidence in the method of locating the critical point by 
monitoring the quantity (s 2 ), even when the system is 
subjected to periodic boundary conditions [lij . 

Exactly the same simulation technique is applied to 
the 3DLG and the results are shown in Fig. [21 The 
finite-size analysis of (s 2 ) in order to extract a is here 
done on the basis of system sizes L =8, 12, 16, 20, 24, 
28, 32, 36 and 40 and up to 2.0- 10 5 MC sweeps are used. 
As in the 2£>CG we see that the system exhibits two dis- 
tinct regimes, one insulating regime consisting of tightly 
bound dipoles and one metallic. There is no scaling of 
(s 2 ) with L below T ~ 0.32 with the system size. Note 
that the change in scaling of the mean square dipole mo- 
ment occurs at a significantly lower temperature than in 
the 2DCG. This is to be expected, since there is more 
configurational entropy available in 3D than in 2D. It is 
worthwhile comparing this result with the one obtained 
in Ref. where the coupling 47r 2 /i corresponds to the 
temperature T here. There the critical value t c = \2tt 2 
was obtained, corresponding to T c = 1/3 in our case and 
agreeing well with our numerical result. 

The main result is that in the 3Z3LG and the 2DCG, a 
low-temperature regime exists where positive and nega- 
tive charges are bound in tight dipole pairs. This regime 
is separated from a high-temperature regime where at 
least a finite fraction of charges are free. The results 
obtained in three dimensions have the same features as 
those found in 2D. The scaling exponent a(T) for (s 2 ) 
is zero in the low-temperature phase and positive in the 
high-temperature phase. Such a change in a(T) can- 
not be analytical, and therefore the two scaling regimes 
must be separated from each other via a phase transition. 
Since (s 2 ) is a measure of the polarizability of the system, 
we conclude that the above demonstrates a non-analytic 
change in the polarizability of the system, i.e., a non- 
analytic change in the dielectric function of the system 
as a function of temperature. Hence, a system of point 
charges with overall charge neutrality interacting with a 
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dard renormalization group (RG) calculation, where in 
a perturbative treatment only analytic singularities are 
generated. Therefore, in a standard RG analysis an ana- 
lytic correction ~ k 2 to the dispersion can be generated. 
On the other hand, since |fc| 3 is non-analytic, it cannot 
be renormalizcd within a standard RG analysis. It is 
conceivable that this is one of the reasons why a recent 
RG analysis of this theory did not find evidence for a 
dielectric phase 0. 

We emphasize that although it is known that a metal- 
insulator transition occurs in 2D via a KT transition 
|l3j| . our numerics by themselves do not demonstrate 
this. To establish the KT nature of the transition on 
purely numerical grounds requires convincing numerical 
evidence that there is a universal jump in the inverse di- 
electric constant at the transition. The main point of 
the present simulations is that they settle the difficult 
question of whether a low-temperature dielectric regime 
exists at all in a 3DLG, which is a system not subject 
to standard electrostatics, and for which the usual the- 
orems on screening of charges by dipoles do not apply. 
The answer is in the affirmative, and we have applied 
our method also to the 2D case for comparison and as 
a benchmark on the correctness of method of monitor- 
ing (s 2 ) to establish the existence of two scaling regimes 
separated by a phase transition. 
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FIG. 2: Results from simulations of the 3DLG. Up to 2.0 ■ 10 5 
MC sweeps were used at each temperature, a) (s 2 ) vs T for 
some of the simulated system sizes L =8, 12, 16, 20, 24, 28, 
32, 36 and 40. Error bars are smaller than the symbols used, 
b) a(T) vs T found from fitting the data of (s 2 ) at different 
T to AL a where A is a constant. A selection of such fits 
are shown in c). It is evident that a low-temperature regime 
where (s 2 ) is independent of L exists in the 3D log-gas in the 
same that that it exists in the 2D log-gas. 



bare logarithmic interaction undergoes a metal-insulator 
transition in both 2D and 31?. 

The 3.DLG can be shown to be equivalent to an anoma- 
lous sine-Gordon theory in d = 3 where the usual k 2 
dispersion is replaced by a |fc| 3 one yj, which is non- 
analytic. This leads to technical difficulties in a stan- 
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